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£f) ' In the framework of a SUSY SO(10) model a phase is generated 

. spontaneously for the B — L breaking VEV. Fitting this phase to the 

' observed CP violating K,B decays all other CP breaking effects are 
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explicitly calculated and found to be in the right range and sign for 
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CP violation is directly observed only in the decays of the K and B mesons. 
The present experimental results are consistent at the moment with the 
standard model (SM). I.e. CP breaking is induced by a phase in the Cabibbo, 
Kobayashi, Maskawa (CKM) mixing matrix of the quarks. 
Extensions of the SM using right-handed (RH) neutrinos, that account for 
the neutrino oscillations, involve in general phases which allow for CP viola- 
tion in the leptonic sector also. This CP breaking is difficult to observe but 
may be detected as soon as neutrino factories are available. The observation 
of neutrino-less double beta decays may be also an indication for Majorana 
phases in the neutrino sector |2j. 

Spontaneous generation of baryon asymmetry in the universe (BAU) needs 
CP violation It is clear now that it requires also extension of the 

SM, while baryon asymmetry in the universe (BAU) a la Fukugita and 
Yanagida [3] due to leptogenesis [H] is the most popular and promising the- 
ory for the BAU. 

Where is the CP breaking coming from? 

CP breaking can be induced via phases in the Yukawa coupling, in the in- 
teractions of the LH and RH gauge bosons and in the VEVs. Phases in 
the spontaneously generated VEVs lead naturally to violation of CP. This 
spontaneous breaking can also help to solve the strong CP problem [Jj jSj. 
The spontaneous violation of CP was already suggested long ago by T. D. Lee UJ. 
In the framework of 50(10) GUT spontaneous breaking was first discussed 
by Harvey, Reiss and Ramond ^Uj. Recently, Bento and Branco added 
to the SM a heavy Higgs scalar with a B — L violating VEV to generate 
spontaneous CP violation. 

In general, the known CP violation in the hadronic sector is not related to 
the leptonic one. Even the CP breaking needed for leptogenesis is usually 
independent of that in the leptonic sector. Hence, CP violation in the lep- 
tonic sector is in general not predictable. Predictability can be gained only 
in terms of a specific model. There are quite a few models relating CP vio- 
lation in the neutrino sector to leptogenesis ^2] but no conventional SUSY 
GUT which connects the leptogenesis to the observed violation in the K and 
B decays is presently known. 

I would like to suggest in this paper that the one and only origin for CP 
violation is a spontaneous breaking at high energies. A phase in the B — L 
breaking VEV can induce all manifestations of CP violation. This phase 
can be fixed by the observed breaking in the K and B decays and the other 
CP violations are then predicted. In particular, we will show explicitly that 
within a SUSY SO (10) model the amount and sign of leptogenesis are as is 
needed to have the right BAU. 
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Let me first show how a phase can be spontaneously generated in the SU(5) 
singlet component of a scalar 126 representation of SO(10). It was already 
pointed out by Harvey, Ramond and Reiss |10j that there is a natural way 
to break CP spontaneously at high energies. This is due to the fact that 
(126) 4 is 50(10) invariant. is the Higgs representation used to break 

down B — L . Its SU(5) singlet component gives also masses to the heavy 
RH neutrinos. The corresponding large VEV induces also small VEVs in 
the components of 3> j^S that transform like 2^ under the SM ^3] which 
play a role in the light fermion mass matrices. 

Assume that all the parameters in the SO(W) invariant Lagrangian are real. 
If the three fermionic families are in \l/i6's, only 3>io, ^T26 an< ^ ^120 can 
contribute to the mass terms: 

16 x 16 = (10© 126) 5 0(12O)as (1) 

Suppose we have chosen global symmetries that dictate a (super-)potential 
of the form 1 jllj 



F(A 1 ,A 2 ,....) = Fo+[...+Al(*lo)|][(*l 2 6)|+(*l26)|]+A 2 [(*126)|+(*l26)|] 

(2) 

and that those are the only phase dependent terms after the spontaneous 
breaking 2 . If the 577(5) singlet component of $126, *T26 ac Q u i re a VEV as 
well as the right component of <&io: 

< *io >= ^= < *126 >= ^=e ia (3) 

The phase dependent part of the potential can be then written as 

V(v,T,a) = Acos(2a) + Scos(4q). (4) 

For B positive and \A\ > 4B the absolute minimum of the potential is 
obtained with 

a = -arccos(— ). (5) 

This spontaneous generation of a phase in the large VEV T, will generate 
also phases in the induced small VEVs which give mass to the light fermions. 
Those will lead to CP violation in the quark and lepton sectors. The value 
of the spontaneously generated phase a depends on arbitrary parameters in 



1 Note that 10 is a real representation. 

2 For a detailed discussion of possible scalar potentials see Ref. The [(3?i2e)s + 

(*T2e)sl P ar t serves also to break the continuous global symmetries avoiding massless 
Nambu- Goldston bosons. 
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the Higgs potential. Its actual value can be however fixed by the require- 
ment that the phases of the induced light VEVs will give the observed CP 
violation in the K,B decays. All other manifestations of CP violation will be 
then explicitly predicted. In particular the amount and sign of leptogenesis 
are then predicted in models where M° %rac is known. 

Let me now explicitly calculate the amount of leptogenesis in a SUSY 50(10) 
model where a phase is generated spontaneously in the B — L breaking VEV. 
The model was developed in a series of papers m] ^S]. It was originally 
aimed to find explicitly the mixing angles which are hidden in the SM, like 
RH rotations. Those allow to calculate explicitly e.g. the proton decay 
branching ratios as well as all mass matrices and in particular the Dirac 
neutrino mass matrix and the RH neutrino mass matrix which are needed 
for the calculation of the leptogenesis. We will use here the mass matrices 
given in Ref. ^Sj. This is a renormalizable SUSY 50(10) model i.e. B — L is 
broken via 3>j26 + <I?i26 while gives mass to the RH neutrinos (without 
using non- renormalizable contributions). The origin of CP breaking in the 
model is a phase in the SU(5) singlet component of one ^j26- A global 
horizontal symmetry U(1)f dictates the asymmetric Fritzsch texture |16| for 
the fermionic mass matrices and the possible VEVs in the different Higgs 
representations. By fitting the free parameters to the observed masses and 
CKM matrix a set of non-linear equations is obtained. These equations have 
five solutions which obey all the restrictions, i.e. five sets of explicit mass 
matrices. The Dirac neutrino mass matrices have the texture: 



~jljf Dirac 






A 





B 





c 





D 


E 



(6) 



They are given explicitly in Appendix I. 

The RH neutrino mass matrices have the following form in our model: 



I a \ 

a \m r . (7) 

y o o -b j 

Where the real a, b > 0. The corresponding eigenmasses are given in Table 1. 



Table 1. The masses of the RH neutrinos for the five solutions in 10 13 GeV. 



Solution 10 13 GeV 


1 


2 


3 


4 


5 


Mi = M 2 


5.2 


9.1 


16 


18 


12 


M 3 


8 x 5.2 


7 x 9.1 


3 x 16 


3 x 18 


2 x 12 



3 



What is Leptogenesis? 

Out of equilibrium CP violating decays of RH neutrinos, Ni, produce excess 
of the lepton number 5L ^ . This will induce baryon asymmetry through 
B + L conserving sphaleron processes @] [3] [H] . 
The amount of CP violation in these decays is: 



Knowing the details of CP violation in the leptonic sector as well as the RH 
mixing angles 3 , one is able to calculate explicitly the BAU via leptogenesis. 
This is the main test of the model. 

Let us denote the Dirac neutrino mass matrix M® zrac in the basis where 
M VR is real diagonal with positive eigenvalues: Md- In this basis €{ can be 
expressed as follows 



T(Nj -> Lj + g) - rCjVjt -> Ljt + $t) 
T{Ni - U + $) + r^t _> L .t + $t) 



1 



^Jm[(MiM D )?.]/(M?/M l 2 ) 



where 



/(x) = + -) + 



x- 1 



2 



and u = 174 x sin/3 GeV 4 . 



is given in eq.(7) and its eigenmasses in Table 1. 



It is diagonalized by a matrix U 



U T M VR U = diag(M 1 ,M 2 , M 3 ) = M 3 diag(j±, ^ 
U = OP where P = e~ ii/2)a diag(i, 1, i) 



1) 



and 




In this basis, in terms of eq. (6) 



3 Note, that M is diagonalized using the RH mixing matrix. 
4 tan/3 = 10 is used in the model |15| 
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M ] D M D 



f 1/2(|A| 2 + |5| 2 + |L>| 2 ) 
-i/2{\A\ 2 - \B\ 2 + \D\ 2 ) 
K l/V2(5Ct - DEfl) 



i/2(\A\ 2 - \B\ 2 + \D\ 2 ) 
1/2(\A\ 2 + \B\ 2 + \D\ 2 ) 



l/y/2(B*C - D^E) 
i/y/2{B^C + D^E) 
\C\ 2 + \E\ 2 



This gives the following general results 

Im((M ] D M D ) l2 (M ] D M D ) l2 ) = Im({M ] D M D ) 2l (M ] D M D ) 2l ) = (8) 

(MjjM^n = (M ] D M D ) 22 

Due to the degeneracy of Ni, N 2 , the decay of both contributes to ej. How- 
ever, eq. (8) avoids the possible singularity in f(x). Hence, 
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-(Im[(Mj,M D ) 2 3 ] + Im[(MlM D ) 2 23 ])f(M 2 /M 2 ). 



8ttv 2 (M' d M d )ii 

The BAU is given then (in the minimal supersymmetric SM) as 

Y B 



-1/3% 

9* 



where g* = 228.75 and ds-L is the dilution factor due to inverse decay 
washout effects and lepton number violating scattering. It must be ob- 
tained by solving the corresponding Boltzmann equation. There are different 
approximate solutions in the literature. The frequently used approximate 
solution 1 7 is good only for 



K 



fhiMp 



fhi (eV) 



1.7 x 87ru 2 V<f 1-08 x 10- 3 (eV) 



> 1 



where 



. In our model however, 



K 10" 



_ (A^M g )ll 
'"1 — Afi 

Buchmiiller, Di Bari and Pliimacher |Hj studied recently in detail both cases 
K > 1 and K < 1. They found that for K < 1 one must take into account 
thermal corrections ue to the gauge bosons and the top quark. Hence, cIb-l 
depends on "initial conditions" and they found that for K ~ 10 -2 5 



10" 4 > 



< io- 



Hirsch and King 18 give empirical approximate solutions for the case 
K <C 1 . The solution corresponding to our model is 

Log w (d B -L) = 0.8 x Logi (rhi eV) + 1.7 + 0.05 x Lo^Mi/lO 10 GeV). 
5 See Figure 9 in their paper where ds-L is called Kf. 
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I will use this expression to have a definite prediction. The results for the 
five solutions are given in Table 2. 

Table 2. The CP asymmetry e^, the dilution factor ds-L and the Baryon 
asymmetry Yb for the five solutions. 



Solution 




d-B—L 




1 


-6.5 x 10" 7 


0.0064 


6.1 x 10~ 12 


2 


-6.6 x 10" 5 


0.0074 


7.1 x 10- 10 


3 


-7.4 x 10" 5 


0.0088 


9.5 x 10~ 10 


4 


-1.3 x 10" 6 


0.009 


1.7 x lO" 11 


5 


-5.6 x 10" 5 


0.06 


4.9 x 10~ 10 



This must be compared with the experimental results: 
BOOMerANG and DASI QH 

0.4 x 10~ 10 <Y B > 1.0 x 10~ 10 

WMAP and Sloan Digital Sky Survey [H] 

Y B = (6.3 ±0.3) x 10" 10 

Hence: 

• Solutions 1 and 3 are probably excluded. The other solutions are consis- 
tent with the experimental observation, especially if the uncertainty in ds-L 
is taken into account. 

• All solution have the right sign. This is the main prediction of the model 
in view of the uncertainty in ds-L- 

I must emphasize that there is no ambiguity in the prediction of the sign 
because of the following reasons: 

a) The sign of Mi must be positive because ej is calculated in terms of Mb, 
which is the neutrino Dirac mass matrix in the basis where the RH neutrino 
mass matrix (7) is diagonal, real and positive. 

b) The parameters and especially the phases of M°"' ac (6) are fixed without 
ambiguity for each one of the above solutions, although one cannot write 
explicitly their dependence on a. As was mentioned before, the entries to 
the mass matrices are solutions of non- linear equations in which the induced 
components of (with the phase a) are involved. The physical value of 
a is then fixed by requiring that Jjariskog ~ 10~ 5 to be a ~ 0.003 6 . 



6 In a recent paper Frampton, Glashow and Yanagida in ref. 1121 presented a model 
where the sign of the BAU can be related to the CP violation in neutrino oscillation 
experiments. In our model both CP violation in the neutrino oscillation as well as the 
sign of the BAU are predicted in terms of CP violation in the quark sector. 
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To complete the predictions of the model let me use the complex lepton 
mixing matrix Upmns of Ref. ^S] (see Appendix II) to give the amount of 
CP violation in the neutrino oscillation 

Jleptons = ^(^11^22^12^21) 

and the value of < m ee > 

3 

< m ee >= ^2(U e i) 2 mi 
i=l 

relevant for the neutrino-less double-beta decay (3f3§ v 7 . See Table 3. 

Table 3. The CP violation invariant for the leptonic sector Ji ep tons an d the 
effective neutrino mass for the neutrino-less double-beta decay for the five 
solutions. 



Solution 


1 


2 


3 


4 


5 


Jleptons 


0.0092 


0.000059 


9.8xl0" 6 


7.8xl0" 6 


6.6xl0" 6 


< m ee > 


0.0031 


0.005 


0.0068 


0.0056 


0.0029 



CONCLUSIONS 

I presented in the paper the following observations: 

CP is naturally broken spontaneously at high energies in 5O(10) GUTs. 

A phase is generated in a VEV and not in the Yukawa couplings, as it is 

usually done. This can be used as the only origin CP violation. 

In the framework of a SUSY 5O(10) model that uses this idea, fitting to the 

observed CP violation, as it is reflected in the CKM matrix, fixes uniquely 

the CP breaking in the leptonic sector without free parameters. An explicit 

calculation of leptogenesis in this model gives solutions consistent with the 

range and sign of the observed BAU 8 . 

Our model applies the conventional see-saw mechanism |22| . it is possible 
however, to use a similar program for the type II see-saw |23| as well |24| . 
The large value of the RH neutrino mass can be incompatible with the 
gravitino problem if SUSY is broken in the framework of mSUGRA. Possible 

7 mi in our solutions is of O(10 -3 eV). 

8 "A common origin for all CP violations" was suggested recently also by Branco, Parada 
and Rebelo (3TJ They use a non-SUSY SM extended by adding scalar Higgs, leptons and 
exotic vector-like quarks. The complex phase is generated spontaneously in the VEV of 
the heavy singlet scalar meson. The connection with the low energy CP violation in the 
hadronic sector is obtained only via mixing with the exotic quarks. They give also no 
explicit value for the leptogenesis. 
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solutions are discussed in the literature. E.g. Ibe, Kitano, Murayama and 
Yanagida presented very recently a nice solution based on anomaly 
mediated SUSY breaking. 
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APPENDIX I 

The Dirac neutrino mass matrices for the five solutions (for tan/? = 10) in 
GeV. 



Solution GeV 


1 


2 


3 


4 


5 


Re(M i f irac )i2 


17.486 


26.953 


-41.320 


-41.320 


-28.274 


Im(M Dirac )i9 

J -- L - L - L \ x V } iz 


0.0394 


0.0607 


0.0929 


-0.0929 


-0.06356 


Re(M^ iroc ) 2 i 


17.654 


27.120 


-41.218 


-41.218 


-28.172 


Im(M^ rac ) 2 i 


0.0394 


0.0607 


0.0929 


-0.0929 


-0.06356 


(M i f--) 23 


-113.425 


-142.425 


116.073 


82.073 


102.073 


(M i f--) 32 


-14.700 


14.302 


10.695 


44.695 


24.695 


Re(M i f- oc ) 3 3 


-127.913 


-176.670 


146.103 


146.103 


78.715 


Im(M i f i ™ c ) 33 


-0.3152 


-0.4249 


0.2788 


0.2788 


0.1271 
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APENDIX II 

The leptonic mixing matrix for the different solutions. 



Solution 


1 


2 


3 


4 


5 


Re([/pM7V,s)ii 


-0.8583 


0.8136 


0.7465 


0.8579 


0.8740 


lm(UpMNs)ll 


0.000004 


0.00034 


-0.000001 


-0.000001 


0.000001 


Re ( Upm N s ) 1 2 


-0.5104 


-0.5778 


-0.6589 


-0.5059 


-0.4806 


lm(UpMNs)l2 


-0.000007 


0.000007 


-0.00027 


-0.00021 


-0.0002 


Re([7 pmns)i3 


-0.0526 


-0.0644 


0.0927 


0.0897 


0.0717 


Im([/pMAfs)i3 


0.000002 


0.00026 


0.00042 


0.00004 


0.00003 


Re ( Upm ns ) 2 1 


-0.3496 


-0.4869 


-0.4653 


-0.3754 


-0.2492 


Im(?7pM7Vs)21 


0.00191 


0.00190 


0.00212 


0.0017 


0.00088 


Re(J7pMAfs)22 


u.uou / 


-U.UiUO 


-u.uiu t 


-U. I ou^ 


-U.tJU 1 u 


Im([/pMAf5)22 


-0.0030 


0.0029 


0.00260 


0.0031 


0.00018 


Re(f/pMAfs)23 


-0.6682 


-0.6185 


-0.6350 


-0.5628 


-0.7829 


Im([/pMAfs)23 


0.0031 


0.00285 


0.0029 


0.0026 


0.0028 


Re(f/pMAfs)31 


-0.3756 


-0.3176 


-0.4756 


-0.3508 


-0.4172 


lm(UpMNs)31 


0.00082 


0.00085 


0.00216 


0.0009 


0.0011 


Re(UpMNs)32 


0.5552 


-0.6168 


-0.4309 


-0.4492 


-0.6664 


lloa(UpMNs)32 


-0.00121 


0.00127 


0.0009 


0.00097 


0.0014 


Re(UpMNs)33 


0.7421 


0.7832 


0.7669 


0.82168 


0.6179 


lm(UpMNs)33 


-0.00163 


-0.00204 


-0.0020 


-0.00213 


-0.0016 
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